o 
o 

(N 



o 



a 



Preprint typeset in JHEP style - HYPER VERSION MCTP-04-38 

INT-PUB 04-17 
hep-th/0406264 

Coupling constant dependence of the shear 
viscosity in AA = 4 supersymmetric Yang-Mills 
theory 



Alex Buchel'^'' 

'^Perimeter Institute for Theoretical Physics 

Waterloo, Ontario N2J 2W9, Canada 

^Department of Applied Mathematics, University of Western Ontario 

London, Ontario N6A 5B7, Canada 



> , , 

Tij" ' Em^ail: abuchel@perimeterinstitute.ca 

\o\ 

CN . James T. Liu 

^^ . Michigan Center for Theoretical Physics, Randall Laboratory of Physics 

■^ , The University of Michigan 



Ann Arbor, MI 48109-1120, USA 



Email: |j imliuSumich . edu 



^- Andrei O. Starinets 

,J^ ■ Institute for Nuclear Theory, University of Washington 

Seattle, WA 98195, USA 



KA , Email: [starinaSphys . Washington . edu 



Abstract: Gauge theory - gravity duality predicts that the shear viscosity of A/" = 4 
supersymmetric SU{Nc) Yang-Mills plasma at temperature T in the limit of large 
Nc and large 't Hooft coupling gy^^c is independent of the coupling and equals to 
7rN^T^/8. In this paper, we compute the leading correction to the shear viscosity in 
inverse powers of 't Hooft coupling using the a'-corrected low-energy effective action of 
type IIB string theory. We also find the correction to the ratio of shear viscosity to the 
volume entropy density (equal to l/47r in the limit of infinite coupling). The correction 
to l/47r scales as (gYj^jNc)'^^'^ with a positive coefficient. 
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1. Introduction 

Transport coefficients such as viscosity, diffusion constants, thermal and electric conduc- 
tivities are the key ingredients in describing the hydrodynamic regime of any medium 
[]T|. These coefficients are usually obtained from experiment rather than computed 
from first principles of an underlying microscopic theory, because the study of realis- 
tic strongly interacting systems often remains beyond the reach of currently available 
theoretical methods. 

For finite-temperature quantum field theories, and thermal gauge theories in par- 
ticular, computations based on the Boltzmann equation in the regime of weak coupling 
g <^1 show 0, 1^, §, 1^ that the shear viscosity behaves as 

riogi/^^ 

where A^^ is the number of colors and C is a large numerical coefficient. 
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Since the entropy scales as 5 ~ A*"^ T^ V3, the ratio of shear viscosity to the volume 
entropy density s = S/V3, behaves as 

"^ "^ >1 (1.2) 



r-^ 



s g^logl/g"^ 

in the regime of weak coupling. On the other hand, hydrodynamic models used to 
describe elliptic flows observed in recent heavy ion collision experiments at RHIC seem 
to favor small values of the ratio r]/s 0, 0, §. This is not a contradiction, however, 
since the experimental data were obtained for the range of energies where the coupling 
remains relatively large, whereas the result ( |1.2| ) is valid for small coupling. 

It is therefore desirable to obtain results for the viscosity-entropy ratio in the 
regime of intermediate and strong coupling. Lattice simulations cannot address the 
issue of real-time dynamics directly facing (among other things) a formidable problem 
of analytic continuation. (For indirect approaches, see [Q, |[TD|, [|n|], ]r2|.) 

In the absence of more conventional methods, the AdS/CFT (or gauge theory /gravity) 
duality conjecture in string theory [|I^, 0, |T^, [T^ emerged as a source of insights into 



the non-perturbative regime of thermal gauge theories. The best studied example of the 
duality relates a four-dimensional finite-temperature Af = 4 SU{Nc) supersymmetric 
Yang-Mills (SYM) theory in the limit of large Nc and large 't Hooft coupling gy^^c 
to the supergravity background corresponding to a stack of N^. near-extremal black 
three-branes. 

Using the AdS/CFT conjecture, one is able to predict the behavior of the entropy 
of A/" = 4 SYM in the regime of strong coupling [|1^, [|18|. In the large Nc limit, the 

entropy is given by 

27r2 
S = —NlV,T'f{glMNc) , (1.3) 

where the function figyM-^c) interpolates (presumably smoothly) between 1 at zero 
coupling and 3/4 at infinite coupling. The strong coupling expansion for / was obtained 
by Gubser, Klebanov and Tseytlin^ [IIH 



figvMNc) = I + ^cmg^YMNc)-'/' + ■■■■ (1-4) 



At weak coupling, / can be found by using finite-temperature perturbation theory |T^ 



3 + V2, 



fig'vMN.) = 1 - ^tAmN. + '-^{gluN^f' + ■■■ (1.5) 



"'^Normalization convention for gvM adopted in the current version of this paper differs from the 
one used in ||lS| . Sec Eq. (p]9) and footnote H. 
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Our goal in the present paper is to obtain the analogue of the strong coupling expansion 
( |1.4]) for the shear viscosity and for the ratio rj/ s. In the limit of infinite N^. and infinite 
't Hooft coupling, the shear viscosity of A/" = 4 SYM plasma was found to be ||20| 

TT 



7] 



-NlT\ (1.6) 



Taking into account the leading term in (|1.4| ), one may notice that [^ 

l = l^^O{llNl)+0{{gl,,N,r'l'). (1.7) 



The shear viscosity (|1.6|) was obtained by a number of interconnected methods [Q, [21 



p2| , p3| all of which relied on the validity of the conjectured gauge/gravity correspon- 
dence. In that sense, the result (|1.6|) does not enjoy the status of a firmly established 
field-theoretical calculation and should be viewed as a prediction of a gauge/gravity 
duality for a particular supersymmetric theory^. 

It was further found in |2^ that the ratio r]/ s remains equal to l/47r for other 
gauge/gravity duals where rj and s alone behave very differently compared to the basic 
example of A/" = 4 SYM. The same feature holds for theories dual to M-branes 



The calculations in |^ were done independently by using the traditional AdS/CFT 
strategy, and also by deriving a generic formula for rj/ s by a method reminiscent of 
the old "membrane paradigm" approach^. It was later shown in ||28| that the type II 
supergravity equations of motion guarantee that the generic value for r]/ s is universally 
equal to l/47r. Whenever the near-horizon supergravity geometry dual to a finite 
temperature gauge theory allows for an extension to asymptotically fiat space-time"^. 



the universality of r]/ s can be related [^] to the universality of low energy absorption 



cross sections for black holes observed in 30 



If the conjectured universality is indeed true, it is most probably associated with 
the regime of strong coupling in various theories. For weakly coupled theories, the ratio 
7]/ s is typically very large, even at large N^. It was conjectured in [^ that 1/4-7? is the 
lowest possible value for rj/ s for all substances in nature. In view of this conjecture, 
finding the first nontrivial term in the strong coupling expansion for rj/ s is especially 
interesting. 

Our paper is organized as follows. In Section ^ we outline the method of computing 
the strong coupling expansion for shear viscosity of A/" = 4 SYM from supergravity. In 

^Over the years, AdS/CFT has survived numerous checks at zero temperature jT^]. For a check at 
finite temperature, see [ p2[ . 

•^From the AdS/CFT point of view, this method presumably computes the lowest (hydrodynamic) 
quasinormal frequency of a given gravitational background. The corresponding link is not firmly 
established at the moment, see pfl], p7| . 

*Such an extension is not known for explicit examples discussed in p8[. 
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Section ^ we consider gravitational perturbations of the a'- corrected black brane metric 
and compute the boundary term of the corresponding on-shell action. In Section ^ we 
use the Minkowski AdS/CFT prescription to find the retarded thermal correlator of a 
stress-energy tensor in the hydrodynamic approximation and compute the correction 
to the shear viscosity. Finally, in appendix ^ we collect the coefficients of the effective 
action which are too cumbersome to appear in the main text. 

2. Shear viscosity of AA = 4 SYM from supergravity 

Shear viscosity is the transport coefficient that appears in the hydrodynamic consti- 
tutive relation for the spatial components of the stress-energy tensor (see |^, |3ll] for 
more hydrodynamics preliminaries). 

Among the methods developed in |20|, [21 1, |23] for computing the shear viscosity 
from supergravity, the most straightforward and the least technically complicated one 



20i is based on the Kubo relation 



r/=lim— / dtd'xe'-'{[T,y{x)T,ym). 



(2.1) 



The Kubo formula is a particular case of the celebrated fluctuation-dissipation theorem 
which relates fluctuations of a medium to the response of the medium to an external 
action. The formula expresses transport coefficients of a slightly non-equilibrium system 
in terms of the real-time correlation functions computed in an equilibrium thermal 
ensemble. Eq. (|2.1|) can be written in the form 



V 



lim 



G 



xy,xy 



{Uj,0) - G^y,^y{ 



^,0) 



^0 2u:i 
where the retarded Green's function is defined as 

-Wd^xe-^«-^^(t)([T^,(x),T,p(0)]) 



G^u,xpi^^Q) 



(2.2) 



(2.3) 



andG^{uJ,q) = {G^{LJ,q))*. 

To compute the hydrodynamic limit of the retarded correlation function ( |2.3| ) from 
supergravity, one can use either the Minkowski AdS/CFT prescription [Q, ||3^, or the 
method based on calculating the absorption cross-section of gravitons by non-extremal 



three-branes [^ |2^. In the absence of a' corrections, to the lowest order in cu/T <^ 1, 
q/T ^ 1 the retarded correlation function is given by ||2ll 



G 



R 

xy,xy 



2rp2 



[^,(1) 



16 



{i2TiTuj + q^ + 



(2.4) 
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From Eq. (|2.2| ) one finds the shear viscosity 

n = iNJT' + O [l/Nl) + O {{gUN.)-'") 



(2.5) 



Corrections to the shear viscosity in powers of the inverse 't Hooft couphng correspond 
on the string theory side of AdS/CFT duahty to a' corrections to classical general 
relativity. The precise relation between the expansion parameters in supergravity and 
Yang- Mills theory is given by^ [16| 



a 
1? 



\/9\m^c 



(2.6) 



Since the leading order stringy correction to type IIB supergravity is proportional to 



a , on the gauge theory side we expect the coupling constant correction to shear 
viscosity to scale as {gyM^c)'^^'^- At the same time, finding 1/N^ corrections requires 
quantum gravity calculations in the black three-brane background. 



In this paper, we use the prescription of ||2^ to compute the retarded correlator 
^xy,xy from the a'-corrected black three-brane metric found in [0. The coupling 
constant correction to viscosity then follows from the Kubo formula 



3. Perturbations of the a'-corrected near-extremal three-brane 
background 

We start with the tree level type IIB low-energy effective action in ten dimensions 
taking into account the leading order string corrections MB, 15^ 



-^^10 J 



R-\{dct>?-^{F,f + ---+ 7e-i<^iy + --- 



1 



(3.1) 



7=gC(3)(a')S (3.2) 

where W may be chosen in an appropriate scheme to be proportional to the fourth 
power of the Weyl tensor |T 

W 



^ ^pmnq^h '"" rsk 



1 



"T 2 ^pqmn^h ^ rsk ■ 



(3.3) 



^ The normalization convention L'^ /a' ^ = AirgsNc = gyM^c IM used in the current version of this 
paper differs from the norniahzation L"^ /a' ^ = 2gYM-^c used in the earher versions and in some other 
works by a factor of 2. In general, this normalization depends on the normalization of the generators 
of the gauge group: g^j^j = 2TTgs/c for tr (TT'') = cS"''. See e.g. Section 13.4 of M. 
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It is helpful to rewrite W in the form |^| 



W = Bijki (25* -^ — _B'J j ^ Bijki = C"^ijnC"'ikm- (3.4) 

In (p.l|) , ellipses stand for other fields not essential for the present analysis. As usual, 
it is assumed that the self-duality condition on F5 is imposed after the equations of 
motion are derived. The form of the action ( p.l|) and subtleties associated with the 
self-duality condition on the five- form are discussed in [|I^, [^. We emphasize that 
the possibility of additional corrections of order {a'Y associated with the five-form is 
not excluded. However, indirect arguments presented in |T^, [^ suggest that these 



corrections, if present, will not affect the near-horizon geometry of the black three-brane 
solution. We proceed with our calculation while ignoring this potential problem. 

3.1 The a'-corrected black three-brane background and its S^ reduction 

In the absence of a' corrections, the metric describing the non-extremal black three- 
brane solution trivially factorizes into the five-sphere part and the "AdS" part. How- 
ever, this is no longer the case for the a'- corrected solutions, as was first emphasized 



by Pawelczyk and Theisen |^ . Thus in considering both the a'- corrected black brane 
solution and its fluctuations, the correct procedure is to work in ten dimensions and 
then perform the S^ reduction to obtain the five- dimensional asymptotically AdS ge- 
ometry. An additional subtlety comes from dealing with the self-dual five-form F5, 
which identically squares to zero by virtue of ten-dimensional self-duality. 

The ten-dimensional metric can be taken in the form of a standard ansatz used for 
the S^ reduction of type IIB supergravity (see e.g. [p8| ) 



dsjo = e-^"(")^5mn(a;)rfx"*rfx" + e^'^'^'^dnl , (3.5) 

where the five-dimensional asymptotically AdS metric g^mni^) has a general form con- 
sistent with the symmetries of the problem^ 

dsl = ^e^(") (-/e'*(")rft2 ^ ^^2 ^ ^y2 ^ ^^2N ^ ^^Ku) (3,6) 

u Au^f 

Here f{u) = l — u"^, tq is the parameter of non-extremality of the black brane geometry, 
and we set the "AdS radius" L to one. When 7 = in ( p.l|) , the metric (|3.5| ) with 
u = 0, a = 0, b = 0, c = is the standard solution of type IIB low-energy equations 
of motion describing the near-horizon limit of non-extremal three-branes. Corrections 



^In a dual thermal field theory, Poincare invariance is broken, but the system is assumed to be 
isotropic. 
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to that solution can be found by solving equations of motion perturbatively in 7. To 
leading order in 7, functions a, b, c, v were found in [1^, [^ 



a(M) = -15 7(5^2 + 5m^-3m^), (3.7a) 

h{u) = 157 (Sm^ + 5m^ - V^u^) , (3.7b) 

c(n)=0, (3.7c) 

u{u) = ^u\l + u') . (3.7d) 



The Hawking temperature corresponding to the metric (|3.5| ) is 



T = — (1 + 157). (3i 



7r 



The dilaton also receives a' corrections |T^. Since to leading order in 7 they do not 



mix with the gravitational perturbations which we are considering in this paper, we do 
not write them explicitly. 

Substituting the ansatz (|3.5|) into the action (|3.1|) we find 

.3 



71-3 f 



R^+2Qe-"-^''-\{d<pf-^{duf-^e-^''+-ie-i^-'-^''W{u,u,-i) 



(3.9) 

where the ten-dimensional gravitational constant kiq = vSvrG is related to the number 
Nc of coincident branes and the AdS radius L by Kiq = 2t\'^^L'^ jN^- 
For small z/ ~ 7 we have 

e~f ''iy(M, z/, 7) = ISOu^ + \%^^u^v{u) + 0{y^) . (3.10) 

(Note that the expression for W{u^ u, 7) is not merely the five- dimensional W computed 
with the metric ( p.6| ).) 

As mentioned above, there is a subtlety in reducing the term involving the five-form 
field strength. The correct procedure is to assume that F5 has components only along 
the five-sphere, and double that contribution in the effective ten-dimensional action. 
Another way of saying this is that the term —8e~^'^ was added to ( ^.9] ) by hand to 
ensure that, upon variation, the correct reduced equations of motion in five dimensions 
are reproduced. 

3.2 Metric perturbations 

The coupling between the boundary value of the graviton and the stress-energy tensor of 
a gauge theory is given by hy{x)T^/2. According to the AdS/CFT prescription, in order 
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to compute the retarded thermal two-point function of the components of the stress- 
energy tensor entering the Kubo formula ( p.2|) , we should add a small bulk perturbation 
hxy{u,x) to the metric (p3|), and compute the on-shell action as a functional of its 
boundary value h 



xy\ 



Simple symmetry arguments [^ show that for a perturbation of this type and 
a metric of the form ( |3.5| ) all other components of a generic perturbation h^^, can 
be consistently set to zero. (We have checked explicitly that to linear order in the 
perturbation, the equation for hxy decouples from all other equations describing a 
generic perturbation of the background (|3.5|) .) 

It will be convenient^ to introduce a field ip{u,x) 



(piu,x) 



2 '^xyy^j *^J 



r, 



(3.11) 







(note that h^ 
sition 



gXXj^ 



xy 



e 3 ip{u^x) 



f{u,x) 



(p{u, x) + 0(u^)) and use the Fourier decompo- 



(2^)' 



-iii)t+ik-x ( \ 



(3.12) 



Since in ( p.2|) we need the correlator at vanishing spatial momentum, it will be sufficient 
to consider perturbations which depend on the radial coordinate and time only. 

3.3 The effective action 

Substituting the perturbed a'-corrected black brane metric into the ten-dimensional 



action (|3.1| ) and performing the five-sphere reduction keeping in mind the subtleties 
mentioned in Section |3?T|, we obtain the effective action for the field (pk{u). To quadratic 
order in ipk{u) and linear order in 7 this action is given by 



where S^^^ is independent oiipkiu), and 



(3.13) 



5(2) 



N^ 



57r" 






du 



A iflip-k + B ip'kV'_k + C ip'kV-k 



+ D v9fc^_fe + E iflip'Lk + F iflip'^ 



(3.14) 



The coefficients A, B, C, D, E, F are even functions of the momentum. They are given 
explicitly in appendix ^ 

''In the absence of a'-corrections the effective action of (p{u, x) is that of a minimally coupled scalar 
in the AdS black hole background. 



Variation of S^"^^ leads to 



6TT^ 



(2vr) 



du [EOM] 5ip_k + {Bi5^-k + B25^'_^) 



J 



where the coefficients of the boundary term are given by 

Br = -{A^k)' + 25^1 + C^k - 2{Eviy + F^l - {F^', 

B2 = Aipk + Fip'^ + 2Eipl, 
and EOM denotes the left hand side of the Euler-Lagrange equation 



Aipl+Cip'^ + 2Dipk-^, {2B^i + dp, + Fv'L) + ^ {A^k + 2E^l + F^',) 



d_ 
du 



(3.15) 



(3.16) 
(3.17) 



0. (3.18) 



In order to have a well-defined variational principle, one has to add a generalized 
Gibbons-Hawking boundary term to the action ( p.l4| ). Variation of this additional 
term is supposed to cancel the contribution B25ip'_^ in ( p.l5| ) without leading to any 
new boundary terms involving ffist and higher derivatives of Lfk- In the absence of 
higher- derivative terms in the action, this is achieved by adding the standard Gibbons- 
Hawking term ||3^ proportional to the trace of the extrinsic curvature of the boundary, 
i.e. by adding the boundary term —AipkV'_^. to the action ( |3.14|) . Similarly, one can 
cancel the second term in E2 by adding the boundary term —Fip'j^ip'_i^/2. However, the 
last term in (|3.17| ) poses a problem: adding a term —2E{p'l{p'_i^ to cancel it would upon 
variation produce second derivatives of Sip^k- This problem is not unexpected because 
we have a theory whose equation of motion is a fourth order differential equation. In 
our case, however, the difficulty can be avoided since we intend to treat the higher- 
derivative term as a perturbation, and we only solve the equation of motion (|3.18 



perturbatively to the first non-trivial order in 7. This means that effectively we are 
still solving a second order equation. Let us write Eq. (|3.18|) in the form 



^l + PiV'k + Po^k = 0(7) 



(3.19) 



where all 7-dependent terms are exiled to the right. Now, if we add to the action (p.l4| ) 
the generalized Gibbons-Hawking term 



k: 



-Aipk(p'_ 



-vW-k + E (pi<^fe + 2pQ^k) <^- 



k 1 



(3.20) 



and note that the coefficient E is proportional to 7, then upon variation we find that 
the boundary term involving 5ip'_f^ vanishes on shell up to terms of order 7^. Thus, to 
linear order in 7, the variational problem is now well-posed. 
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The bulk action ( p.l4|) can be rewritten in the form 



5(2) 






where 



A' C 

F , , F' , 



(3.22) 



and EOM denotes the left hand side of (p.l8|) . Adding to this the generalized Gibbons- 
Hawking term ( p.20|) , we observe that the complete on shell action reduces to a bound- 
ary term 



(2^ 



^k 



(3.23) 



where 



^k 



87r2 



{B - A) i^'kip-k +2^^- A') ifkf^k - ^VfcV-fe + E^W-k 



F' 



^ + ^\j^j 



to^ 



Wk^-k - -7rvW-k - E :^vW-k + "^E—^if'j^if.k 

Z UJ UJ 



(3.24) 



and t» = uj/2ro. 

3.4 The solution for h 



xy 



We now turn to finding the solution to the equation of motion (|3.18|) . Explicitly, 
Eq. (|31^ ) reads 



- 7 — ^ (-600 + 6193^2- 25552 m^ + 22327m^ + 592 M^tt)^) 0'^'(r) 

- 7 — ^ (-600 + 819 M^ - 3440 u^ + 7105 u^ + 2368 u^ tt>^)0fc(r) 



u 



7 



8/ 



100 u - 550 u^ + 1000 M^ - 450 u^ + 600 1»^ - 719 u^ w^ + 3968 m^ vo 



4^2 



2665 M^ tt)^ + 4 w^(74t»4 - 25) 



<Pk{r). 



(3.25) 
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where we put all 7-dependent terms to the right. Eq. (|3.25| ) can be solved perturba- 
tively in 7 by writing 

^u{u) = vt\u)+^^^^\u), (3.26) 

where ip]^ \u) is the solution to the minimally coupled massless scalar equation (i.e. 
Eq. (lOSi ) with 7 = 0) found in ^, 



^t\n) 



[l-u] 2 ( 1 - — In— — + 0(t»^ 



(3.27) 



The full solution obeying the incoming wave boundary condition at the horizon u = 1 

(3.28) 



2T| and normalized to one at the boundary m = is 

ipk{u) = {l-uy~G{u), 
where the function G{u) is regular at -u = 1, and is given explicitly by 



G{u) = 1-7- 



16 






16 



log(l + u) 



+ 7- 



iwu 



(43 u^ + 135 u^ + 195) + 0(7^ w 



(3.29) 



4. Coupling constant correction to shear viscosity 

Having found the solution for a gravitational perturbation, we can compute the corre- 
lation function Gxy,xy{^,(l) by applying the Minkowski AdS/CFT prescription ^Q\ 



qR 



xy,xy 



[uj,q) = lim 2jFfc , 



u^O 



(4.1) 



where J-'k is the boundary term (|3.24|) . (It is assumed that all contact and momentum- 
independent terms in (|4.1| ) should be discarded.) Computing the limit, to the lowest 
order in w we find 



U^O 4:71'' 



1 + 757 - ^ - ivo{l + 1957) + 0(i»2 



M" 



(4.2) 



On a technical side, we note that the asymptotic behavior of the coefficients ( [A.l| )- 
(|A.6|) and the regularity of the solution (|3.2(j|) at u — > imply that only the first two 



terms in ( p. 241 ) contribute to the limit ( |4.2| ), with the nontrivial contribution coming 
only from the term {B — A) ip'j,ip_k. Moreover, contributions to the coefficients A and 
B coming from the Weyl term W vanish in the limit. The W term therefore plays the 
role of a "grey cardinal" influencing the result indirectly through the correction to the 
black brane metric and the equation of motion. 
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Thus the retarded correlator at zero spatial momentum and to the leading order 
in ui is 

G%,.y{uj, 0) = -t -^ (1 + 1957) = -t —^ (1 + I5O7) , (4.3) 

where we used the expression (|3.8|) for the Hawking temperature of the a'-corrected 
metric. Applying the Kubo formula (|2.2|) , we immediately find 

r^ = '^N^T^ (1 + I5O7) . (4.4) 



Combining ( p.6|) and ( p.2|) , we obtain 



,8 



V = In^T' (1 + ^ C(3) (^L/iVc)-^/^ + ■■■). (4.5) 



Using the result of Gubser, Klebanov and Tseytlin ( |1.4|) for the entropy density we find 
the strong coupling expansion for the ratio r]/s 

- = ^ (1 + 1357) . (4.6) 

s Att 

Thus for large 't Hooft coupling Qy^^c ^ 1 the correction to the ratio of shear viscosity 
to the entropy density in A/" = 4 supersymmetric Yang-Mills theory is positive, 

^ = ^(l + ^C(3)(^?-MiVe)-^/^ + ■■■). (4.7) 



Formulas ( ^4.5| ) and ( [4.7| ) are the main result of this paper. Thus we have shown that 
the conjecture of [|23[ that t]/s > l/Air remains valid to the next order in the 't Hooft 
couphng expansion. 
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A. Coefficients of the effective action 

The coefficients of the five- dimensional effective action ( |3.14|) are given by 



A 



B 



C 



2u 
8u 



8f{u) + 7^2 (-600 + 25^2 + 1760m^ - 1185u^ - Uu^w'^) 



(A.l) 



24/(n) + 7^2 (-1800 + I79u^ + 5424m^ - 3387u^ - 7Q8u^w'^) 



(A.2) 



4^2/ 



8/(m)(3 + u^) + 7^2 600 - 825u^ - 16945u^ + 34005u^ - 16835m* 



+ 104M^t»2 + 872M^t»2 



(A.3) 



ArpA 



D 



8u^P 



Idfiuf + 8uf{u)w^ + -fu^ 250n + 1570n^ + 2220n^ - 1920m^ 



+ 600t»2 + 25uW - 2400M*t»2 + 10Q7u''w^ + 8M'^(37t»* - 265) 

^ = 7377^^TV/(M)^ 

F = 727r^TV/(M)(ll-37M2). 



(A.4) 

(A.5) 
(A.6) 
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